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ABSTRACT 

Finding and characterizing the population of active galactic nuclei (AGNs) 
that produces the X-ray background (XRB) is necessary to connect the history 
of accretion to observations of galaxy evolution at longer wavelengths. The year 
2012 will see the deployment of the first hard X-ray imaging telescope that, 
through deep extragalactic surveys, will be able to measure the AGN population 
at the energies where the XRB peaks (~ 20-30 keV). Here, we present predictions 
of AGN number counts in three hard X-ray bandpasses: 6-10 keV, 10-30 keV and 
30-60 keV. Separate predictions are presented for the number counts of Compton 
thick AGNs, the most heavily obscured active galaxies. The number counts are 
calculated for five different models of the XRB that differ in the assumed hard 
X-ray luminosity function, the evolution of the Compton thick AGNs, and the 
underlying AGN spectral model. The majority of the hard X-ray number counts 
will be Compton thin AGNs, but there is a > 10 x increase in the Compton thick 
number counts from the 6-10 keV to the 10-30 keV band. The Compton thick 
population show enough variation that a hard X-ray number counts measurement 
will constrain the models. The computed number counts are used to consider 
various survey strategies for the NuSTAR mission, assuming a total exposure 
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time of 6.2 Ms. We find that multiple surveys will allow a measurement of 
Compton thick evolution. The predictions presented here should be useful for all 
future imaging hard X-ray missions. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert - 
surveys — X-rays: diffuse background — X-rays: galaxies 



Introduction 



The vas t majority of the growth of supermassive black holes (SMBHs) is d riven by ac 



cretion (e.g.. Marconi et all 120041 : Merloni fc Heinzll2008l : iDi Matteo et al.ll2008f ). Therefore 



a complete census of accret ing SMBHs throughout cosmic time is necessary to quantify the 
efficie ncy of accretion (e.g.. lSoltanlll982l ; lYu fc Tremaindl2002t IShankar et al.ll2004t lYu fc Lu 
20081 ). whic h, in turn, elucidates the connection between black ho l e growth and galaxy evo- 
lution (e.g.. lHaiman et alJliool l2007t Eopkins k Hernquistlbood : IShankar et al.ll2009h . All 
efficiently accreting SMBHs are intrinsically luminous X-ray sources, while the underlying 
host galaxy is, in general, a much fainter X-ray emitter. The cosmic X-ray Background 
(XRB) is then naturally interpreted as being the integral emission of all accreting SMBHs in 
the universe, and the hard slope of the XRB spectrum indicates that most of the AGNs are 
obscured behind substantia l columns of gas and dust in the host galax y (e.g. , ISetti & Woltjer 



1989; Comastri et al. 1995; Ueda et al 



20031 : iGilli et all 1200 ll . 120071 ). This expectation was 



spectacularly confirmed by the numerous deep surveys performed by Chandra and XMM- 
Newton in the 2-10 keV band that resolved > 80% of the XRB into individual objects 



(e.g. 

tions of these sour ces show that a 
redshifts of z ~ 1 ( Alexander et al. 



Brandt fc Hasingerl 120051 ; IWorsley et al.ll2005l ). Multiwavelength follow-up observa 



20021 ; iBarger et ajj 12005] ; iBrusa et all 120101 ) 



arge fraction of them are obscured, type 2, AGNs at 



2001 



Tozzi et alJbOOll IBarger et al.1 12002 IStern et al. 



However, summing the emission from all the AGNs detected by Chandra and other mis- 
sions working below 1 keV cannot account for the observed peak of the XRB at ~ 30 keV 
(e.g. . IGilli et al.l 120071 ) . The missing AGNs are likely dominated by Compton thick sources, 



those behind obscuring column densities Nn > o"t ~ 10 cm . The combination of pho- 
toelectric absorption and Compton scattering reduces the observed flux of these AGNs to 
such a degree that, if they lie at any reasonable cosmol ogical distance, they are invisible 
even in the deepest Chandra or XMM-Newton surveys (IBrandt fc Hasingerl 120051 ) . There 
are currently three strategies that are being employed to discover and identify Compton 
thick AGNs. The first is to look for the 'waste heat' being re-radiated by the absorber 
in the mid-infrared. Several groups have uncovered a population of Compton thick candi- 
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2006; Daddietal. 2007; Alexander et al 



dates at z > 1 using th i s idea (IStern et al J 12005c lAlonso-Herrero et al.ll2006c iPolletta et al. 



2008 : 



Fiore et al.l 



2008 



2009), but the constraints 



are compl icated by the necessity to separate out contributions f r om obscured star forma - 
tion (e.g.. lMushotzkyll2004l : bonlev et alll2008[ feckart et alfeoiot beorgakakis et al.l 120 101 ) . 
The second method to detect Compton thick AGNs is to search at hard X-ray energies 
(> 20 keV) where the sources are less affected by obscuration (although still suffer from 
Compton scattering losses). The hard X-ray detectors on board BeppoSAX, Swift, and IN- 
TEGRAL have been able to reveal local (i.e., z ~ 0) Compton thick AGNs thr ough this 



method (e.g.. iBassani et al.lll999t IVignali fe Comastril |2002| ; iMalizia et al.ll2009al ). Recent 



surveys by the latter two missions have measured the fraction of all AGNs that are Compton 
thick at z ~ to be ~ 10-2 0%, depending on the luminosity range considered and the 



details of the source selection (Sazonov et al 



Treister et al. 2009: Winter et al. 



2009 



2007 



Burlon et al. 



Tueller et al. 2008; Malizia et al. 2009b 



201l|). Finally, Compton thick AGNs 



may also be identified i n pointed X-ray observations by searching for objects with specific 
optical char acteristics (jGilli et al.l 120111), or those that have very large Fe Kg equivalent 



width s (e.g., Levenson et al. 20061 ; LaMassa et al. 2009 ; Comastri et al. 2011 ; Feruglio et al. 



201ll ). However, these observations require very long exposure times and the translation 
from Fe Ka line measurement to column density dep ends on model-dependen t details such 



as the geometry and metallicity of the absorber (e.g. iMurphy fe Yaqoobll2009l ). 



The hard X-ray instruments flown in previous missions have had fairly poor sensitivity 
when compared to detectors at lower energies. Therefore, there is very little known about the 
Compton thick AGN population at z > 0. Intriguingly, if one takes the infrared estimates of 
Compton thick densities at high re dshift at face value, then there must be strong ev olution 
of the Compton thick population (IDraper fe Ballantynd |2010| ; iTreister et al.ll2010l ). This 
evolution might indicate that 'Compton thickness' is an unique phase of AGN evolution and is 



tied t o specific Eddington ratios or merger events (IDraper fe Ballantyndl2010l ; ITreister et al. 



20101 ) . Therefore, discovering the missing Compton thick AGNs beyond the local Universe 
and tracing their evolution is vital to understanding a host of problems related to black hole 
fueling and galaxy evolution. 

In 2012 the Nuclear Spectroscopic Telescope Arra^ (NuSTAR; Harrison et al. 2010 ) 
will be launched as part of NASA's Small Explorer Program. This mission will be the first 
ever focusing hard X-ray telescope covering an energy range of 6-80 keV. It will therefore 
be imaging the sky at the energies where the XRB peaks and the Compton thick sources 
are most visible. One of the primary objectives of the baseline two-year NuSTAR mission 



1 http : //www. mis tar . caltech. edu/ 
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is to perform extragalactic surveys to describe the AGN population and XRB emission at 
these energies. Other hard X-ray imaging missions, such as Astro-H (Takajiashi et al. 2010 ; 



launching in 2013/2014) and the proposed New Hard X-ray Mission (jTagliaferri et al.ll2010l ). 
will follow later this decade. The data produced by these missions will provide crucial tests 
of models of AGN and galaxy evolution. 

Thus, in light of the upcoming launch of NuSTAR, there is a need for models of AGN 
evolution, and, specifically, Compton thick evolution, that can be used to compare against 
future hard X-ray data. In addition, model predictions can also be used to assist the planning 
of the hard X-ray extragalactic surveys to ensure that a reasonable number of Compton thick 
AGNs are detected at several different redshifts. To be the most effective, several models, 
spanning a range of parameters, should be available to compare against the survey data 
with the goal that many of them will be eliminated. As a first step, this paper predicts 
the hard X-ray number counts for three different measurements of the AGN hard X-ray 
luminosity function with each requiring a different Compton thick fraction. In addition, 
we also predict counts for a model where the Compton thick fraction is a function of AGN 
Eddington ratio, and therefore varies with redshift and AGN luminosity. Finally, we also 
present a conservative model that predicts the smallest fraction of Compton thick AGNs 
consistent with the available data. Our results are specifically designed for comparison with 
the upcoming hard X-ray surveys, and are the first to show the effects of the different model 
parameters on the AGN and Compton thick counts. 

The next section describes the various XRB models used to calculate the Compton 
thick and total AGN number counts. The counts are presented in Section [3j Sect. H]then 
applies these results to the design of deep hard X-ray surveys, with a specific example for 
NuSTAR. Our conclusions for the NuSTAR mission and for other future hard X-ray missions 
are summarized in Section This paper assumes the following cosmological parameters: 
H = 70 km s" 1 Mpc" 1 , Q A = 0.7, and Q m = 0.3. 



2. Calculation of Predicted Number Counts 

The integrated AGN number counts, N(> S), in a specific energy band are 
N(>S) = — / ; ' 77- wo ' ^ d\ogL x dz, 

H oJz min ymax(logL x (mm),logL x (S NH )) d\ogL X (1 + Z) 1 (Q m (1 + Z) 6 + S2 A ) 1/2 

(1) 

where d&(L x , z)/d\ogL x is a given hard X-ray luminosity function (HXLF), di is the lu- 
minosity distance to redshift z, and L x (Sn h ) is the unabsorbed rest-frame 2-10 keV lu- 
minosity that gives an observed frame flux Sn h for a source at z averaged over the Nu 
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distribution for that Lx and z. The integrals are evaluated from z m - m = to z meLX = 5, 
and \og(L x (min)/erg s -1 ) = 41.5 to log(Lx(max)/erg s -1 ) = 48. The integrated counts are 
computed in three energy banda^: 6-10 keV, 10-30 keV and 30-60 keV. As described below, 
the calculation of the expected number counts relies on a number of quantities which are 
poorly constrained or subject to significant uncertainties. 



2.1. Unabsorbed AGN Spectrum 



The conversion between luminosity and flux requires a spectral model for the average 
AGN population at a given Lx and z. The assumed unabsorbed rest-frame AGN spectrum 
is comprised of a power-law with photon index T that varies from source to source (e.g., 
Winter et al.l 120091 ) . an exponential cutoff at E cut = 250 keV, and a reflection component. 
The distribution of the power-law cutoff energies is not observationally constrain ed and will 



be an important measurement for the upcoming hard X-ray missions. However, iGilli et al. 



(120071 ) find that including a range of cut-off energies in their X-ray background synthesis 
model has a negligible effect on the results; therefore, we do not consider this parameter to 
be a significant source of uncertainty in the predictions. 

The reflection component hardens the spectrum above > 10 keV and produces an Fe Ka 
line. A reflection spectrum may arise from eithe r the accretion disk or from cold, dense 
gas a t large distances from the black hole (e.g., iRoss fe Fabianl 120051 ; iMurphy fe Yaqoob 
20091 ). Thus, the strength of the reflection component in an AGN spectrum may vary 
signifi cantly from source to source, althoug h it is known that, on average, it cannot be very 
large jGandhi et al.ll2007l : lBallantvne![2^1of l. In addition, there is observational evidence that 
the strength of the reflection spect r um from the di stant reprocessor is anti-correlated with 
the luminosity (IBianchi et al.l 120071 ; IShu et al.l 120101 ) . 



2.2. Ah Distribution and Compton Thick AGNs 

The basic spectral model will suffer some obscuration that depends on the fraction of ob- 
scured AGNs at a specific L x and z, as well as the distribution of obscuring column densities. 
The fraction of absorbed, or type 2, AGNs (f'2) is defined as those objects observed through 
column densities Ah > 10 22 cm -2 . This fraction is observed to decrease with AGN luminosity 



2 To assist comparisons with previous observations and models, Appendix A shows predictions in the 
5-10 keV band. 
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tentatively, increase with z (ILa Franca et al 



e.g., lUeda et al.l 120031 : Barger et al.ll2005l: lLa Franca et al.ll2005l: iTueller et al.l 120081). and 



2005 



Ballantyne et al 



2006 



Treister fc Urry 



20061 ; lHasingerll2008l ). The distribution of column densities that comprise the obscured and 
unobscured populati on has been measured in the local univer s e, or in a redshift and luminos- 



ity integrated sense (IComastri et al.lll995t iRisaliti et al.lll999t lUeda et al.ll2003l ; iDwelly et al. 



20051 ; llbzzi et al.ll2006l ). However, it is possible (perhaps likely) that the distribution of col- 
umn densities will also depend on L x and z. 

Previous deep X-ray surveys by Chandra and XMM-Newton are not sensitive to most 
Compton thick AGNs; therefore the space density and evolution of these objects represents 
the largest uncertainty in the census of AGN. Additional unobserved objects have to be 
added to the HXLF in order to boost the integrated emi ssion to fit the XRB spect rum, 
and/or the measured z = Compton thick space density (IDraper fc Ballantyne] l2009l ). As 
mentioned above, Compton thick AGNs are obscured by columns Ah > 10 24 cm -2 , but once 
Ah reaches a few times 10 24 cm -2 the transmitted continuum becomes severely suppressed 



flMatt et al.l Il999h . Indeed, extreme Compton thick AGNs with A H > 10 25 cm" 2 will be, 
for all practical purposes, invisible even in the deepest hard X-ray surveys. Consequently, 
such objects cannot provide a significant contribution to the XRB, and this paper defines 
Compton thick AGNs as those with columns 24 < log Ah < 25. 

The simplest method for including Compton thick AGNs into a XRB model is as a 
straightforward extension of the Compton thin population that evolve with other obscured 
AGNs (i.e., along with f 2 ). This is quantified by the Compton thick fraction, / CT , defined 
as the fraction of all type 2 AGNs that are Compton thick. For example, /ct — 0.5 means 
that there are as many Compton thick AGNs as Compton thin type 2s. 



2.3. X-ray Luminosity Functions 

In recent years there have been several measurements of the HXLF that have signifi- 



HXLF ( 


Ueda et al. 


2003 




Bareer et al. 


2005 




La Franca et al. 


Ebrero et al. 


20091: 


Yencho et al. 


2009; 


Aird et al. 


2010 


). For 



Ah distribution, these differences in the HXLF can result in factors o f 2-3 variations in the 
number density of Compton thick AGNs (IDraper fc Ballantyne! 120091 ) . 
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2.4. Model Definitions 

Given all the uncertainties associated with the calculation of number counts, we con- 
sider five separate models to determine the impact of our ignorance of many of the above 
parameters. The most interesting physical parameters that will likely be constrained by the 
hard X-ray number counts are the HXLF and Compton thick fraction. Therefore, we con- 
sider three XRB models with a fixed Compton thick fraction that differ only in t he assumed 



HXLF . The fourth XRB model uses one of the test HXLFs (specifically, the lUeda et al. 



(120031 ) HXLF), but assumes a more complicated form for the Compton thick evolution, 
where the Compton thick fraction depends on the Eddington ratio. The fifth and final 
model proy i des a lower-limit to the number of Compton thick AGNs. This model uses the 



Ueda et al.l (J20Q3J) HXLF and a fixed Compton thick fraction, but assumes a brighter AGN 
spectrum and the minimum allowable Compton thick density. Thus, this model predicts the 
most conservative view of the Compton thick population that is consistent with the available 
dated 



2.4.1. HXLF Variable Models 



In the following models, the reflection spectrum is calculated using the 'reflion' model 
(IRoss fe Fabian! 120051 ) assuming a solar Fe abundance. This reflection spectrum is averaged 
over all viewing angles, and is added to the power law such that the equivalent width o f 
the narrow Fe Ka line agrees with the observed X-ray Baldwin effect (IBianchi et al.l 120071 ). 
The corresponding reflection fraction drops from ~ 1 at low luminosities to close to zero 
at high luminosities. To account for the observed distribution of photon indices, both the 
power-law and the reflection s pectrum are Gaus sian averaged about T = 1.9 with r m i n = 1.5, 
r max = 2.3 and a r = 0.2 (e.g., billi et al.ll2007h . 



The fraction of obscured AGNs, f 2 , is assumed t o be dependent on both redshift and L y , 
such that /2 oc (1 + z) a (\og Lx)~ b where a = 0.40 ( Ballantyne et al. 2006 ; Treister fe Urry 
20061 ) and b = 4.7. The normalization of f'2 is set by fixing the type 2 to type 1 ratio at 4 : 1 



interested readers can perform a similar test by computing the number counts for the lGilli et al.l (|2007l ) 
XRB synthesis model using their webtool at http://www.bo.astro.it/~gilli/counts.html. Note that 
the Compt on thick fra c tion o f this model is similar to the HXLF variable models used here. As Appendix A 
shows, the iGilli et al.l (|2007j) predictions are within the range of those studied in this paper, and do not 
change any of our conclusions. 

20081 ). then the 10-30 keV Compton thick counts are 



If a = 0.62 and continues out to z = 2 (Hasinger 



increased by < 13% at 1 < z < 2, but are reduced by < 9% at z < 0.5. 
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for z = and log(Lx/erg s 1 ) = 41.5, as determined from observa tions (jMaiolino fc Rieke 



19951 ) and XRB modeling flBallantvne et all 120061 : iGilli et all 120071 ). The redshift evolution 
is halted at z = 1 corresponding to the flat tening evolution of the cosmic star formation 
rate density (e.g., iHopkins fc Beacoml 120061 ). The unabsorbed type 1 sources are evenly 
distributed over the columns logiV# = 20,20.5,21, and 21.5, while Compton thin type 2 
sources are evenly distributed over the columns log Nh = 22, 22.5, 23, and 23.5. 

The Compton thick fraction, /ct, is determined by requiring the Compton thick AGN 
space density at z ~ with log(Lx/erg s _1 ) > 43 to be w 3.6 x 10~ 6 Mpc~ 3 , in agreement 
with the a bsorption corrected local space density measurement from Swift /BAT and IN- 
TEGRAL jTreister et al.lbood ; iRigbv et al.lbood: IPraper fc Ballant vnel boioh. This value is 



very similar to the density predicted by the IGilli et al.l (l2007f) model (e.g., I Alexander et al. 
20081 ). which, judging from Fig. 3 in the paper by lTreister et al.l (120091 ). suggests that ~ 15% 
of z ~ AGNs in this luminosity range are Compton thick. The Compton thick AGNs are 
evenly distributed over the columns logiVff = 24, 24.5, and 25, consistent with the limits 
from local observations (ITreister et al.ll2009Q. and t h e Com pton thick absorption profiles are 
taken from the Monte-Carlo models of iMatt et al.l (119991 ) . A scattered reflection spectrum 
with a luminosity of 2% of the intrinsic luminosity is used in place of the transmitted ab- 
sorption spectr um for the log N h = 25 spectral model, and is added to the logiVn = 24.5 
spectrum (e.g.. IGilli et al.l 120071 ). 



The three HXLFs used in these models (in order of increasing Compton thick fraction) 
are the ones published by Ueda et al. (2003; /ct = 0.3, blue lines in all figures), La Franca et 
al. (2 005; / ct = 0-4, green lines) an d Aird et al. (2010; /ct = 0.5, red lines). The lUeda et al. 



(120031 ) and lLa Franca et al.l (120051) HXLFs are best modeled by luminosity dependent den- 
sity evolution (LDDE), while lAird et al.l ( 120101 ) find that luminosity and density evolution 
(LADE) is the best description of their data. 



2.4.2. Compton Thick Evolution Model 



In addition to the above three models, we also consider the scenario where the Compton 
thick fraction (/ct) evolves sep arately from the less obscur ed AGN. Specifically, we make 
use of the composite model of Draper fc Ballantynd (120101 ) who assumed that /ct wa s a 
function of the AGN Eddington ratio. After consider ing several scenarios and co mparing to 
the high-z, IR-derived Compton thick space densities, Draper &: Ballantynd (120101 ) estimated 
that Compton thick AGNs make up ~ 86% of all AGNs with Eddington ratios greater than 
0.9, ~ 60% of AGNs with Eddington ratios less than 0.01, a n d 0% of AGNs at intermediate 



Eddington ratios. This model makes use of the lUeda et al.l (120031 ) HXLF and recovers the 
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observed z = Compton thick density (IDraper fe Ballantyndl2010l ). All other parameters 
are identical to the previously described models. This model is shown as the black line in 
all figures. 



2.4-3. The Conservative Compton Thick Model 



We define a model, described by lTreister et al.l ( 120091 ). that provides a lower- limit to the 
Compton thick population that is consistent with all available data. This model assumes a 
constant power-law with T = 1.9 and a cutoff energy of 300 keV, atte nuated by photoelectric 



absorp tion and a constant Compton reflection component as given by lMagdziarz &: Zdziarski 
( 19951 ). The reflecting material is assumed to provide a reflection fraction of 1.0, has twice 



the solar iron abundance, and has an average inclination angle to the line of sight of 45 deg. 
The larger iron abundance and constant reflection fraction result in a slighter more luminous 
AGN spectrum above 10 keV when compared to the models described in Sect. l2.4.Tl and l2T4"72"l 



The A% distribution employed here is the one shown in figure 2 of iTreister et al.l (120091 ) 
and was originally derived from a torus geometry corresponding to an AGN unification 
model with an aspe ct ratio matching the local ratio of Type 2 to Type 1 Seyfert galaxies 
( Treister et al.112004 ). However, it is fully c onsistent with the observed distributions in various 



hard X-ray survey s 



Comastri et al.lll995l : iRisaliti et al.lll999l : lUeda et al.ll2003l : iDwellv et al. 



2005 ; Tozzi et al. 20061 ) . In this Ah distribution the majority of the Compton thick AGNs 
have iVfi < 10 24 5 cm -2 , which results in a more luminous average Compton thick spectrum 



than used in the models described above. Since this paper is focused on either the Compton 
thick AGNs, or on the total AGN population, the remaining differences in the two Ah 
distributions have no impact on the results. The luminosity dependence of the obscured 
\<-'y iv.w-j-;,,,, r i^o ar iy from /2 = 1 to between Lx = 10 42 and 3 x 10 46 erg s _1 



(Treister & Urrv 


2005) 


the 


Ueda et al. 


( 


2003|) 



HXLF is assumed. 



The Compton thick fraction of this model is set to be consistent with the observed 
fraction of AGNs at z ~ that are Compton thick, as deduced from hard X-ray sury eys by 



Swift and INTEGRAL (lSazonov et al.ll2007l : iTueller et al.ll2008l : ITreister et a 
corre sponds to a fraction of all z ~ AGNs that are Compton thick of ~ 10% ( 



2009). This 



Treister et al. 



20091 ; Fig. 3), and is likely to be a strict lower-limit due to the observational biases of those 
surveys. The extinction corrected density used to normalize the Compton thic k fraction in 



the m odels described in Sect. 12.4.11 attempts to partially address those biases ( IRigby et al. 



20091 ). Thus, this model presents the most conservative view of the Compton thick popula- 
tion, and is shown as the cyan line in all plots. 
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2.5. Summary of Models 



Predictions of the hard X-ray number counts rely on several parameters or distributions 
that are only partially observationally constrained. Thus, five separate models were defined 
to fully explore the sensitivity of the results to these assumptions. Of the five, three ex- 
plored differences in the measured AGN HXLF (Sect. 12.4. ip . one considered a model where 
the Compton thick fraction was a function of AGN accretion rate (Sect. I2.4.2p . and one 
provides the most conservative view of the Compton thick population (Sect. |2~4~73]) . Figured] 
shows the X-ray background spectrum, as well as the contribution from Compton thick 
AGNs for these five models. All models provide an adequate fit to the X-ray background 
spectrum, which nicely illustrates the degeneracy of synthesis modeling. The redshift evo- 
lution of the log(Lx/erg s -1 ) > 43 Compton thick space density predicted by the models 
is shown in figure |2j The differences in the three different HXLFs are easily ob served with 
the C ompton thick density increasing with the necessary /ct- At z < 1. the lUeda et al. 
(120031 ) HXLF requires the least amount of Compton thick AGNs, and the lAird et al. I (bold ) 



HXLF h as the most. Ab o ve z > 1, t he different evolut ions of the HXLFs become apparent 
and the lLa Franca et al.l (120051 ) and lAird et al.l ( 120101 ) HXLFs predict a similar Compton 
thi ck density. Th e brigh ter spectral model and lower Compton thick fraction employed in 
the lTreister et al.l (120091 ) model results in a low Compton thick density at all z. Finally, the 
evolving Compton thick model predicts a small fraction of high luminosity Compton thick 
AGNs at z < 1 because there are very few high Eddington ratio AGNs to host such obscured 
objects. Thus, there is a significant difference in the predicted densities of Compton thick 
objects at z < 1, implying that deep hard X-ray surveys may be capable to distinguish 
among these models. 



Results 



3.1. Number Counts 



Figures El H] and E] show the predicted N(> S) distributions for multiple redshift ranges 
in the 6-10 keV, 10-30 keV and 30-60 keV bands, respectively. The solid lines plot the 
counts for all AGN, while the dashed lines present the predictions for Compton thick 
source s. The colors d i fferen tiate between t he various HXLFs and Compton thick evolu- 
tions: lLa Franca et al.l (hood) H XLF (green) . lAird et al l 120101) HXLF (red) the composite 
model of Draper fc Ballantynd ( 20101) (black), a n d the lTreister et al.l ( 20091 ) model (cyan). 



For clarity, the model bas ed on the lUeda et al. 



plots. Recall that both the iDraper fc Ballantynd (120101 ) and the iTreister et al.l (120091 ) make 



(2003) HXLF has been omit t ed fro m the 
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use of the lUeda et all fl2003h HXLF 



3.1.1. All AGNs 



Concentrating initially on the counts of all AGNs (solid lines), we see that the numbers 
are v ery similar to those observed in the softe r energy bands by C handra and XMM-Newton 



(e.g., iRosati et al.l l2002r iLehmer et al.l l2005t iMateos et al.l 120081 ) . This is not surprising: 
most AGNs are only moderately obscured and so would appear in both 2-8 keV and 10- 
30 keV surveys. These results imply that most of the sources found in hard X-ray surveys 
will have Chandra and/or XMM-Newton counterparts with similar fluxes, especially for 
AGNs at higher z where the K-correction has moved the absorption out of the Chandra or 
XMM-Newton band (see Sect. Q . 



The lUeda et al.l ( 120031 ) and lLa Franca et al.l (120051 ) HXLF predict very similar numbers 
for the total counts, except for z > 2, where there are very few constraining datapoints. In 
contrast, the lAird et al.l (120101 ) HXLF consistently predicts smaller values for the counts. 



especi ally at fluxes greater than 10 13 erg cm 2 s 1 . This result is because the lAird et al. 



( 120101 ) HXLF predicts fewer high luminosity AGNs at all redshifts. Moreover, the LADE 
model does not require an increase in the space densi ty of low- z, low-luminosity AGNs. The 
predicted number counts from the lAird et al.l (120 lOh HXLF are different enough from the 
other models that surveys will be able to test the viability of this HXLF at different redshift 
ranges. 



3.1.2. Compton Thick AGNs 

Turning now to the Compton thick objects, we see that the number counts are uniformly 
small, comprising only 5-10% of the total AGN counts. Therefore, discovering and measuring 
the evolution of these objects will be a tough observational challenge. However, figures |3H5] 
immediately show the advantage of hard X-ray imaging to find Compton thick AGNs. At 
all redshifts, the number counts of Compton thick source increase by factors of ~ 10-100 
when moving from the 6-10 keV to the 10-30 keV band. This result emphasizes that hard 
X-ray imaging is an important method to identify significant numbers of Compton thick 
candidates. Interestingly, there is little difference in the expected Compton thick number 
counts between the 10-30 keV and 30-60 keV bands. In fact, the counts drop slightly in 
the 30-60 keV band at z > 1 as the K-correction causes the observed frame flux to drop. 
This is a result of the the high-energy cutoff in the AGN power-law and electron recoil in 
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the reflection spectrum, both of which reduce the intensity of the total spectrum at high 
energies. We conclude that, in order to maximize the sensitivity to Compton thick AGNs, 
high energy observations need only to extend out to 30-40 keV (depending on the redshift 
of interest). 

The combination of the different HXLFs and assumptions of the Compton thick evolu- 
tion result in wide range of predictions for the number counts of these AGNs. Figures |3H5] 
indicate that the best discrimination of the models will result from objects at z < 1. In 
particular, the model of iDraper fc Ballantynd ( 120101 ) where Compton thick objects are com- 
prised of a combination of low Eddington ratio objects and very high Eddington ratio AGNs 
(black lines), predicts a significantly different number count distribution than the other four 
models where Compton thick objects are an unchanging fraction of the AGN population, 
and evolve in the same wa y as the less obscured sources. For example, at z < 0.5 the 
Draper fc Ballantynd ( 120101 ) model predicts ~ 2x more Compton thick A GNs with a 10- 



30 ke V flux larger than IP"* 3 erg c m" 2 s" 1 than the lLa Franca et all (120051 ) HXLF (Fig. gj. 



In the IDraper fc Ballantynd (120101 ) model, these high flux, low-z Compton thick AGNs are 
dominated by low accretion rate objects similar to the Circinus galaxy. Pushing to higher 
redshifts moves a sur yey closer to the quasar era where high Eddington ratio accretion is 
more common and the IDraper fc Ballantynd (120101 ) model predicts additional Compton thick 
sources. The factor between the evolving and fixed Eddington ratio models at a depth of 
1Q-13 er g cm -2 g -i j g therefore boosted to ~ 5 in the 0.5 < z < 1 bin. In addition, the 
evolving Compton thick model has a shallower slope at bright fluxes than the fixed Comp- 
ton thick fraction models due t o low luminosity, weakly accreting Compton thick AGNs. 
Thus, the IDraper fc Ballantynd ( 120101 ) model can be easily tested with a Compton thick 
AGN number counts measurement that moves beyond z ~ 0.1 and to a 10-30 keV depth of 
10~ 13 erg cm" 2 s _1 . In this way hard X-ray surveys will be able to measure the evolution of 
Compton thick AGNs and determine their relation, if any, to stages of galaxy evolution. In- 
deed, confirming that the Compton thick AGN fraction is larger at higher redshift will have 
implications for understanding the physics of AGN triggering and the A GN-galaxy connec- 
tion. Finally, as seen by comparing the conservative iTreister et al.l ( 120091 ) model (cyan lines) 
to the other predictions, the number of Compton thick AGNs detected by hard X-ray surveys 
will help test the assumptions made in that model, in particular the shape of the average 
spectral model and the value of the local Compton thick density. 
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3.2. Identification and Optical Detectability of the Hard X-ray Sources 

As seen above, hard X-ray imaging is extremely efficient in detecting AGNs, but ad- 
ditional analysis will be required in order to begin to use these detections to investigate 
problems in galaxy and black hole evolution. This is particularly true for the Compton 
thick AGNs, as they cannot be identified from hard X-ray imaging alone — additional data 
is necessary to determine the level of obscuration for each hard X-ray source. Moreover, 
understanding the nature of the detected hard X-ray sources will require that the host 
galaxies be sufficiently bright so that redshifts and stellar population measurements can 
be made at optical and infrared wavelengths. Both of these problems can be addressed if 
the hard X-ray surveys are performed on fields with significant multiwavelength coverage. 
For example, we have seen that the majority of hard X-ray sources should have Chandra 
or XMM-Newton counterparts provided that at least moderately deep (i.e., a flux limit of 
~ 10 -16 -10~ 15 erg cm -2 s" 1 ) soft X-ray observations exist in the observed field(s). Fig- 
ure E] shows the relationship between the soft and hard X-ray fluxes for the spectral model 
described in Sect. 12.4.11 The left panel plots 0.5-10 keV fluxes for AGNs with (intrinsic, 
unabsorbed) Lx = 10 43 erg s _1 against the 10-30 keV flux for different obscuring column 
densities. As an example of the typical Chandra sensitivity, the horizontal line shows the 
50%-completeness flux limit from the Chandra observations of the All-wavelength Ext ended 



Groth Strip International Survey (AEGIS, eight 200 ks exposures; lLaird et al.l 120091 ) . The 
solid line shows the fluxes for unabsorbed AGNs, and, for all realistic flux limits reached 
by the next generation of hard X-ray surveys (~ 10~ 14 -10~ 13 erg cm" 2 s" 1 in the 10-30 
keV band), these sources will have 0.5-10 keV fluxes greater than the flux limit for mod- 
erate/deep Chandra observations in many survey regions. Thus, ~ 90% of the hard X-ray 
detections will have soft X-ray counterparts that can be combined with the hard X-ray data 
to estimate N^. 

The dashed line plots the fluxes for Compton thick AGNs, and the effects of absorption 
are clearly seen in the predicted 0.5-10 keV fluxes. Although Compton thick AGNs will be 
an order of magnitude brighter in the 10-30 keV band, ones detected in the hard X-rays 
still seem to be at or above the detectability limit of the deepest Chandra surveys (if they 
lie at low z). However, this plot neglects the steep decline in the ACIS-I effective area at 
energies greater than 5 keV. This effect is illustrated in the right panel of Fig. [6] which 
plots the 5-10 keV flux versus the 10-30 keV flux. This figure clearly shows that some 
Compton thick AGNs will be detectable in hard X-ray surveys, but beneath the sensitivity 
limit of deep Chandra surveys. These results again indicate the power of hard X-ray imaging, 
with no decline in effective area above 5 keV, to discover the most heavily obscured AGNs 
between z = and 1. Thus, in addition to those that might be found by combining Chandra 
and hard X-ray detections, Compton thick candidates can also be identified by searching 
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for objects that have no 5-10 keV Chandra counterparts. Analysis using infrared data (e.g., 
Alexander et al.ll2008f ) can then be pursued on these candidates to confirm the identification. 
Finally, there may be an unknown number of heavi ly obscured objec ts that have extremely 



weak scattered components in the soft band (e.g., lUeda et all 120071 ). Such objects would 



also be detected out to moderate redshift by hard X-ray imaging. 

Redshifts are typically obtained from spectroscopic or photometric measurements of 
the host galaxy. Therefore, it is important to consider the expected i?-band magnitudes for 
different 10-30 keV fluxes to obtain a sense of the difficulty of identifying the host galaxy of 
the hard X-ray detected AGNs. This is done by using th e empirical relations bet ween X-ray 
luminosity and absolute i?-band magnitude reported by iLa Franca et al.1 ( 120051 ). measured 
independently for obscured and unobscured AGN: 



logL fl = (0.959 ± 0.025) logL x + (2.2 ± 1.1), 



(2) 



for unobscured AGN, and 



log L R = (0.462 ±0.026) log L x + (23.7 ± 1.1) 



(3) 



for obscured ones, where Lx is the intrinsic luminosity in the rest-frame 2-10 keV band. 
These equations are obviously an approximation, as they do not take into account some of 
the more detailed properties of the host galaxies, like evolving stellar populations or varying 
amounts of extinction^. Such considerations are particularly important for obscured and 
Compton thick AGN, in which the optical light is dominated by the host gala xy. The results 



are in distinguishable from the calculations, based on a dusty torus model, by iTreister et al. 
(120041 ). which assumed the same host galaxy for all AGN. 



For the deepest hard X-ray surveys with an limiting 10-30 keV sensitivity of 10 14 erg cm 2 



most AGNs will have bright optical magnitudes, R < 23 (AB), and, thus for surveys in well 
studied fields, the vast majority of the sources will have detections at other wavelengths and 
measured redshifts. For a very shallow X-ray survey (~ 10~ 13 erg cm -2 s" 1 depth), most 
AGNs will be brighter than R ~ 20 mag (AB), and hence within reach of the Sloan Digital 
Sky Survey. 



5 However, the lLa Franca et al.1 (|2005l) sample has a similar z distribution as will be found in the next 
generation of hard X-ray surveys, so the observed evolution of i?-band luminosity should be similar. 
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4. 



Survey Strategies for Hard X-ray Missions 



4.1. Area vs. Depth 



Deep extragalactic surveys require significant investments of observing time and so need 
to be carefully planned in order to maximize the scientific return. One important consid- 
eration in survey planning is the trade-off between depth and £1X6 cl, clS the luminosities and 
redshifts of sources detected in a very deep and narrow field is often very different from those 
uncovered in a shallower but wider field. The predicted number counts shown in Figs. [3rE3 
contain all the information necessary for survey planning, but require additional manipula- 
tion to see the effects of surveys of differing areas. Therefore, Figs. [7| and [8] present examples 
of a new way of plotting the predicted AGN and Compton thick AGN counts that is specif- 
ically designed for planning hard X-ray surveys. Each figure plots contours of numbers of 
AGNs (upper panel for all AGNs; lower panel for Compton thick AGNs) as a function of 
survey area and 10-30 keV flux. The different colors indicate the predictions for four of 
the five different models considered in this paper (see Sect. I2.4p . The dashed horizontal 



lines in each plot show the area of five regions wit 
Deep Wide-Field Surv ey (9.3 deg 2 ; Murray et ali 



20071 : Elvis et al.ll2009r i. AEGIS (0.5 deg 



i goo d multiwavelength c overage: Bootes 
2005T). COSMOS (2 deg 2 ; Irlasinger et al. 



Davis et al. 



2007 



Laird et ajj|2009j) , the Extended 



Chandra Deep Field South (ECDF-S, 0.25 deg 2 ; iLehmer et aDbood) . and the Great Obser 



vatories Origins Deep Survey (GOODS, 0.089 deg 2 : Alexander et al.ll2003l : iLuo et all 120081 1. 
Finally, b oth panels of figu re [7] also have vertical lines indicating the flux needed (estimated 
using the lUeda et al.ll2003l model) to reach a specific percentage of the integrated 10-30 keV 
X-ray background spectrum. 

These figures contain a large amount of information and should be useful guides for 
planning future har d X-ray sur v eyfl Close examination of the contours shows that, with 



the exception of the lAird et al.l (120101 ) HXLF, all the models predict very similar numbers 
for the total number of AGNs. There is a slight spreading out of the contours at small 
fluxes, illustrating the uncertainty in the slope and evolution of the faint end of the HXLF. 
In contrast, there is a significant difference in the number of Compton thick AGNs predicted 
at a given flux and survey area, indicating that deep surveys will be able to easily distinguish 
between the various models. 

Identifying highly embedded AGNs over a range of redshifts will be one of the principle 
goals of any hard X-ray survey; however, these AGNs will be faint and will require a signifi- 
cant investment of observing time to obtain a sizable number of sources. Fig. |H] shows that an 



'Plots for different energy bands and redshift ranges are available by contacting the authors. 



- 16 - 



area the size of the COSMOS field must be surveyed to a depth of 2-3 x 10~ 14 erg cm -2 s _1 
in the 10-30 keV band to uncover at most 10 Compton thick AGNs at z < 0.5. Similar 
numbers of Compton thick sources would also be found for that depth and sensitivity at 
0.5 < z < 1 and 1 < z < 2. In contrast, a narrow field (such as the GOODS field) would 
have to be observed to a flux limit 10 x fainter to obtain similar numbers of z > 0.5 Compton 
thick AGNs. This survey, however, would not detect as many z < 0.5 objects. 



4.2. Application to NuSTAR 

As mentioned in Sect. [TJ only a portion of NuSTAR's baseline two year mission will be 
dedicated to extragalactic surveys. The limited amount of time that is available requires that 
careful pre-launch planning be performed to optimize the survey strategy. Here we make use 
of the hard X-ray number counts presented above to predict the results of a NuSTAR survey 
of five different regions with significant ancillary multiwavelength data, including deep low- 
energy X-ray observations: Bootes, COSMOS, AEGIS, ECDF-S and GOODS. Each survey 
has a total exposure time of 6.2 Ms, which, assuming 50% efficiency, would take 6 months 
to complete. 

The NuSTAR sensitivities were derived using the "nustar_sens" program provided by the 
NuSTAR instrument team. The program incorporates the NuSTAR effective area, corrected 
for aperture stop, scattering from surface roughness, figure error, detector efficiency, and 
all known attenuation along the optical p ath in the form of ther mal covers and beryllium 



windows. Combined with the GEANT4 (lAgostinelli et al.l 120031 1 simulated internal back- 



ground and the diffuse X-ray background (for a low Earth orbit inclined at 5°), the program 
calculates a sensitivity map for the entire detector given a specific signal to noise, point 
spread function (PSF) and input spectrum. A T = 1.7 spectrum, app ropriate for a mod- 



erately obscured A GN uncorrected for absorption and reflection (e.g., iMalizia et al.l 12003 



Winter et al.ll2009f ). is assumed for all calculations, and we require a 4cx detection. Since the 
code was run for long exposures the signal to noise is in the Gaussian regime. The PSF used 
for the calculations follows a Gaussian+King profile and has an half-power diameter (HPD) 
of 41". Because extraction from a region larger than the HPD will add more background 
than source counts, the HPD is used as the extraction region. 

By definition, surveys require mapping areas greater than the detector field-of-view 
(13' x 13'), so survey sensitivities were calculated for two tiling strategies. The first, called 
'corner shift', is when the survey region is covered by non-overlapping fields of view with 
additional observations at the corners between fields of view. This strategy has the advantage 
that individual exposures can be fairly long, but there is a significant variation in sensitivity 
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over the survey region as NuSTAR, like many other high-energy instruments, has an effective 
area that is a strong function of off-axis angle. The second strategy is called 'half shift', and 
was calculated by shifting the detector sensitivity map over by half a field of view in both 
directions. Each individual pointing is shallower than the corner shift, but the half shift 
provides a more uniform sensitivity map over the survey region. One by-product of the half 
shift survey is that, due to the fact that the pointings are calculated such that they entirely 
cover the given survey size with the best sensitivity, the actual survey size ends up being 
larger than the input, with the extraneous area having lower sensitivity. 

The area versus sensitivity curves for the five surveys are shown in Figure [9] with the 
corner shift surveys plotted as solid lines and the half shift surveys shown as dashed lines. The 
details of each survey — including the number of pointings, average sensitivity, and resolved 
fraction of the X-ray background — are listed in Table [TJ As expected, the corner shift 
surveys probe slightly deeper than the half shift ones, but the latter provide more uniform 
coverage over the survey area, especially in the narrow-fields like GOODS or ECDF-S. The 
numbers of AGNs predicted from both tiling strategies are nearly equal, but the wider area 
provided by the half shift survey increases the AGN yield by ~ 10%. Thus, we conclude 
that half shift tiling should be used for all NuSTAR surveys. 

Table |2] shows the number of AGN detections predicted by each of the five XRB models 
for each of the five half shift NuSTAR surveys. The significant decrease in effective area 
with energy is manifested in the numbers of detectable AGN: the majority of sources will 
be detected in the 6-10 keV band with the 30-60 keV band yielding very few objects and 
basically no Compton thick objects. However, ~ 100 AGNs will be detected in the 10- 
30 keV including (depending on the model) ~ 10-20 Compton thick sources. Moreover, these 
objects will be detectable up to and beyond z ~ 2. These surveys will therefore increase the 
number of AGNs detected in hard X-rays by over an order of magnitude, and will be able 
to estimate the evolution of Compton thick AGNs above z ~ 0. In addition, good spectral 
measurements of absorbing column density, reflection strength and high energy cutoff will be 
determined for a large number of AGNs over a range of luminosity and z. Such measurements 



were only previous ly available for a small number of nearby objects (IBeckmann et al.1 12009 



Molina et al.ll2009f ). Interestingly, because of the higher effective area and the negative K- 
correction, most Compton thick AGNs will also be detected in the 6-10 keV band unless 
the exposures are very shallow (e.g., see the numbers for Bootes field in Table |2]). As the 
NuSTAR effective area peaks at ~ 10 keV, some fraction of the low- z heavily obscured AGNs 
will not be detected by the deep Chandra surveys (see Fig. [6]). 

The choice of NuSTAR surveys must be driven by the science goals: discovering the 
sources that contribute to the peak of the XRB, and determining the evolution of highly 
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obscured sources. The predicted number of Compton thick detections is roughly constant for 
all the surveys (except for the narrowest GOODS-like field), with the redshift distribution 
shifting to higher redshift as the field narrows. This fact allows some flexibility in designing 
a survey; for example, a 6.2 Ms survey with an area ~ 2 deg 2 will detect 15-20 Compton 
thick AGNs with redshifts z < 2. Such a survey will also resolve ~ 30% of the XRB in 
10-30 keV band. The variation in predicted numbers of Compton thick objects implies that 
a simple counting experiment that can push to z ~ 1 may provide important constraints on 
the models. A very deep field such as the ECDF-S will also be of value, as it will probe 
a different part of the ^-luminosity plane, and therefore would increase the statistics of the 
high-;?, highly obscured population, and improve the resolved fraction of the XRB at high 
energies by detecting a larger number of lower luminosity AGNs at high z. 

Finally, recall that the numbers discussed here assume that 6.2 Ms are set aside for ex- 
tragalactic surveys in the baseline NuSTAR two-year mission. However, the sp acecraft has no 



consu mables and is expected to have an orbital lifetime greater than 5 years ( Harrison et al. 



20101 ) . An extended NuSTAR mission will be able to spend several months per year building 
up and extending the surveys performed during the baseline mission. Thus, NuSTAR may 
end its mission with extragalactic surveys that have exposure times 2 to 3 times larger than 
the 6.2 Ms assumed here, with the numbers of detected AGNs correspondingly larger. 



5. Conclusions 



The next few years will see the launch of the first generation of hard X-ray focusing 
telescopes. By focusing X-rays with energies > 10 keV these missions will increase by orders 
of magnitude the sensitivity of observations in this energy range. This capability will be 
vitally important to study the growth of accreting black holes, as much of this accretion 
is expected to be obscured by significant columns of gas and dust. This paper presents 
predictions for AGN number counts in three hard X-ray bands (6-10 keV, 10-30 keV, and 
30-60 keV) to understand and quantify the potential of extragalactic surveys performed by 
these future mission^. The calculations pay particular attention to the expected number 
of Compton thick AGNs, as these heavily embedded objects may be physically connected 
to rapid galaxy growth and are missed by previous observations. We present predictions 
for five different models of the XRB and Compton thick evolution. Specifically, results are 
presented for three different measurements of the AGN HXLF which differ significantly in the 
predicted evolution of high- z low-luminosity AGNs. A fourth model assumes the Compton 



7 The results can easily be extended to other energy and flux ranges by contacting the authors. 
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thick fraction is a function of the AGN Eddington ratio. The final model presents the most 
conservative view of the Compton thick population. These five models predict numbers of 
Compton thick AGNs that differ up to a factor of ten, therefore any measurement of the 
Compton thick number counts beyond z = will highly constrain the fraction and evolution 
of Compton thick AGNs. 

Most hard X-ray detected AGNs will be Chandra or XMM-Newton sources and have host 
galaxies with R < 23; therefore, sources which are well separated from the background should 
be easily identified and have counterparts in multiwavelength ancillary data. Combining the 
hard and soft X-ray data will also allow the Nr estimates that are are necessary to identify 
Compton thick AGNs. Other Compton thick candidates will be identified by following up 
sources that are undetected in the Chandra or XMM-Newton data. In order to make use of 
these data, hard X-ray surveys should be performed in well-observed regions with deep and 
uniform multiwavelength coverage. 

Specific predictions for five extragalactic surveys were performed for the NuSTAR mis- 
sion, scheduled for launch in early 2012. All the surveys assumed a total available exposure 
time of 6.2 Ms, and a half shift tiling strategy is recommended. These deep surveys will yield 
hundreds of AGNs at z < 2, including (depending on the area covered) 10-20 Compton thick 
AGNs over that redshift range. Thus, the NuSTAR deep surveys will allow a detailed explo- 
ration of the nature and evolution of nuclear obscuration in galaxies. The simulations were 
appropriate for the two-year baseline mission of NuSTAR, but the surveys will ultimately 
probe much deeper during the extended NuSTAR mission. 

The predictions presented here show that, assuming equal sensitivities, hard X-ray sur- 
veys can more efficiently detect AGNs and follow their evolution than lower energy X-ray 
observations. It is therefore hoped that future technology development will increase the sen- 
sitivity of hard X-ray instrumentation. A future highly efficient hard X-ray imaging mission, 
in combination with a ground based 30m optical/IR telescope, JWST and ALMA, would be 
extraordinarily powerful in understanding galaxy and black hole growth throughout cosmic 
time. 
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of the National Aeronautics Space Administration under contract NAS8-03060. The authors 
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A. Predicted 5-10 keV Counts 



The hardest energy band that Chandra and XMM-Newton have measured AGN nu mber 
counts is 5-10 keV (e.g jRosati et al.ll2002[ bella Ceca et al.ll2004J : brunner et aDl2008h . For 
completeness, and to help compare against the available data, Figure IA1I plots predictions 
of the total and Compton thick AGN 5-10 keV counts for four of the models discussed 
in the paper, as well as the iGilli et al.l (120071 ) model (magenta lines). The Compton thick 
counts from this last model are bracketed by the conservative model (not shown) and the 
Ueda et al.l ( 120031 ) model at all fluxes. The cyan d ata show the measur ed counts obtained by 
the XMM-Newton survey of the Lockman Hole (IBrunner et al.l 120081 ) . an d the blue points 



plot t he counts from the XMM-Newton Hard Bright Serendipitous Survey (IDella Ceca et al. 
20041 ). 
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Energy (keV) 



Fig. 1. — The solid lines sho w the X-ray backgr ound spectra p redicted by the five AGN evolu- 
tion mo dels conside r ed her e: lUeda et all d2003h HXLF fbluel.lLa Franca et all tood ) HXLF 
(green), lAird et al.l (120101 ) HXLF (red), the iDraper fc Ballantvnd f|2010h composite model 
(black), and the conservative Compton thick model ( Treister et al. 2009^ (cyan). The dashed 
lines plot the Compton thick contribution to the X-ray background. The X-ray background 
data are from the following instruments: blue: ASCA GIS (Kushino et al. 2002); magenta: 
Rossi X-ray Timing Explorer (RXTE; Revnivtsev et al. 2003); green: XMM- Newton (Lumb 
et al. 2002); red: BeppoSAX (Vecchi et al. 1999); yellow: ASCA SIS (Gendreau et al. 1995) 
cyan: XMM- Newton (De Luca & Molendi 2004); grey data: HEAO-1 (Gruber et al. 1999) 
blue data: INTEGRAL (Churazov et al. 2007); red data: Swift /BAT (Ajello et al. 2008) 
black data: Swift/XRT (Moretti et al. 2009); green data: INTEGRAL (Tiirler et al. 2010). 
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Fig. 2. — Predicted space density of Compton thick AGNs with log(L y /erg s 1 ) > 4 3 as a 
function of z. The different models a re distinguished by the line c olor: lUeda et al.l (120031 ) 



HXLF (blue). lLa Franca et al.l fcOOfih HXLF (green), lAird et al.l (1201(1 ) HXLF (red), the 



Draper fc Ballantvnel ( Woh composite model (black), and the conservative Compton thick 



model ( Treister et al. 2009h (cyan ). The triangle denotes the local Compton thick density 



calculated by iTreister et al.l (120090 . while t he circle shows the same density after correcting 



for the flux-luminosity relation reported by iRigby et al.l (120091 ) . 



-28 - 




-14 -13 -12 -14 -13 -12 



log(S) [erg cm 2 s ! ] log(S) [erg cm 2 s 



Fig. 3. — The solid lines plot the integrated number counts of all AGNs in the 6-10 keV 
energy band. The contribution to the counts by Compton thick AGNs are shown as dashed 
lines. The top-left panel plots the counts over the redshift range z = to 5. The other panels 
show the integrated counts over smaller ranges of redshift. The colors diffe rentiate between 



the va rious HXLFs an d/ or Compton thic k evolutions assu med by the model: lLa Franca et al. 
J2005I ) HXLF (green). Eird et al.l (boioh HXLF (red), the braper fc Ballantvne! (boioh com- 
posite model (black), and the conser yative Compton th ick model ( Treister et al.ll2009h (cyan). 



For clarity the model that uses the lUeda et al. 



(120031 ) HXLF has not been plotted. 
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Fig. 4. — As in Figure [3j but for the 10-30 keV energy band. 
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Fig. 5. — As in Figure [31 but for the 30-60 keV energy band. 
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Fig. 6. — Expected fluxes in the 0.5-10 keV band (left) and the 5-10 keV band (right) as a 
function of the 10-30 keV fluxes for AGNs with intrinsic, unabsorbed 2-10 keV luminosities 
of L x = 10 43 erg s _1 . The solid line plots the fluxes for AGNs with N H = io 21 ~ 21 - 5 cm' 2 and 
the dashed line indicate fluxes for AGNs with Nh = io 24-24 - 5 cm -2 . Redshift increases from 
upper right to lower left; the filled circles show the redshifts for each track at z — 0, 1, 2 and 
3 from right to left respectively. The horizontal lines indicate the flux limit of the Chandra 
survey (ei ght 200 ks exposu res) of the AEGIS field to which 50% of the total survey area is 
complete ( jLaird et al.l 120091 ) . 
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Fig. 7. — (Top) Contours of number of AGNs as a function of survey area and 10-30 keV 
sensitivity. Res ults for the model s are plotted using the same color scheme as previous 
figures with the lUeda et al.l (120031 ) HXLF omitted for clarity. These contours include all 
sources in the range < z < 5. The contour levels increase by factors of 10 from 0.1 sources 
(lower-right corner of the figure) to 10,000 sources (upper-left corner of the figure). Only 
the black contours are labeled. The horizontal lines show the areas covered by some well 
known multiwavelength surveys. The vertical lines indicate the fl ux level required t o reach 
a certain percentage of the X-ray background (as judged by the lUeda et al.l 120031 model) 
in the 10-30 keV band. (Bottom) Contours of Compton thick AGNs. The contour levels 
increase by factors of 10 from 0.01 sources (lower-right corner of the figure) to 1000 sources 
(upper-left corner of the figure). 
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Fig. 8. — As in Fig. however these contours include only those sources in the range 
< z < 0.5. 
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Fig. 9. — Area versus 10-30 keV sensitivity for 4<r detections of AGNs in different NuSTAR 
surveys with a total exposure time of 6.2 Ms (about 6 months, assuming 50% efficiency). 
The solid lines are for a corner shift survey, and the dashed lines plot the results for a half 
shift survey. Results are shown for surveying the following regions: Bootes (blue), COSMOS 
(red), AEGIS (black), E-CDFS (green) and GOODS (cyan). 
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Fig. Al. — Predicted total (solid lines) and Compton thick (dashed lines) AGN number 
counts on the 5-10 keV band. The colors denote the diffe rent XRB models as described 
in Fig. [TJ The magenta lines show the predictions of the iGilli et al.l (120071 ) model. The 
cyan crosses are the measured counts from the XMM-Newton survey of the Lockman Hole 
(iBrunner et al.l 12008), and the blue s tars are from the XMM-Newton Hard Bright Serendip- 
itous Survey (IDella Ceca et al.l 12004 ). 
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Table 1. NuSTAR Survey Parameters 



Depth Pointings Avg. Sensitivity (erg cm 2 s 1 ) X-ray Background Fraction 

(ks) 6-10 keV 10-30 keV 30-60 keV 6-10 kcV 10-30 keV 30-60 keV 



Bootes (9.3 deg 2 ) 



Corner Shift 
15.9 392 

Half Shift 

7.4 841 



2.0 x 10~ 14 9.9 x 10~ 14 7.0 x 10 



2.1 x 10" 14 1.1 x 10~ 13 7.2 x 10 



-13 



-13 



0.16 



0.16 



0.14 



0.14 



0.05 



0.05 



COSMOS (2 deg 2 ) 



Corner Shift 

86.4 72 8.6 x 10~ 15 4.3 x 10~ 14 3.1 x 10~ 13 0.28 
Half Shift 

36.8 169 9.4 x 10~ 15 4.8 x 10~ 14 3.3 x 10~ 13 0.28 



0.23 



0.20 



0.07 



0.07 



AEGIS (0.5 deg 2 ) 



Corner Shift 
346 18 

Half Shift 

127 49 



4.4 x 10~ 15 2.3 x 10" 14 1.6 x 10 



5.3 x 10~ 15 2.7 x 10~ 14 1.8x10" 



-13 



13 



0.48 



0.42 



0.35 



0.27 



0.09 



0.09 



ECDF-S (0.25 deg 2 ) 



Corner Shift 
778 8 

Half Shift 

250 25 



3.2 x 10" 15 1.6 x 10" 14 1.2 x 10 



4.0 x 10~ i5 2.0 x 10~ 14 1.3 x 10 



-13 



-13 



0.54 



0.48 



0.35 



0.35 



0.13 



0.13 
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Table 1 — Continued 



Depth Pointings Avg. Sensitivity (erg cm 2 s 1 ) X-ray Background Fraction 

(ks) 6-10 keV 10-30 keV 30-60 keV 6-10 keV 10-30 keV 30-60 keV 



GOODS (0.089 deg 2 ) 


Corner Shift 
3110 


2 


1.8 x 10~ 15 


9.7 x 10~ 15 7.2 x 10~ 14 


0.65 


0.50 


0.15 


Half Shift 
691 


9 


2.9 x 10~ 15 


1.4 x 10~ 14 9.7 x 10~ 14 


0.54 


0.50 


0.13 



Note. — The total exposure time for each survey is 6.2 Ms, constituting a 6 months survey (assuming 
50% efficiency). The average sensitivity is the flux limit of the survey at 50% coverage (i .e., half the 



total a rea). The fraction of the X-ray background at that sensitivity is estimated using the lUeda et al 



(|2003l ) model (see Sect. El 
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Table 2. AGN Detections for Half Shift NuSTAR Surveys 



Model 



Redshift Range 



6-10 keV 



All AGNs 
10-30 keV 



30-60 keV 



Compton Thick AGNs 
6-10 keV 10-30 keV 30-60 keV 



Bootes (9.3 deg 2 ) 



Ueda et al. (2003 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



284 
113 
87 
52 



110 

65 
24 
10 
1 











La Franca et al. (2005 



< z < 5 
< z < 0.5 
0.5 < z < 1 

1< z < 2 

2 < z < 5 



280 
80 
89 
63 
17 



99 
47 
25 
11 
2 



Aird et al. f2010 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



338 
91 
155 
66 
2 



123 
56 
47 
10 




10 

3 
5 
3 




15 

10 

4 

1 













DB10 



< z < 5 
< z < 0.5 
0.5 < z < 1 

1< z < 2 

2 < z < 5 



296 
142 
92 
54 



126 
83 
31 
11 
1 



17 

5 
7 
4 
1 



19 
11 
6 
1 




1 
1 







Treister et al. f2009 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



250 
102 
78 
58 
11 



107 

66 

28 

12 





4 
3 
1 
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Model 



Table 2 — Continued 



Redshift Range 



All AGNs 



Compton Thick AGNs 



6-10 keV 10-30 keV 30-60 keV 6-10 keV 10-30 keV 30-60 keV 



COSMOS (2 deg 2 



Uedaet al. (2003 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



217 
59 
75 
60 
10 



79 
36 
24 
12 
1 



7 
4 
2 
1 












La Franca et al. (2005 



< z < 5 
< z < 0.5 
0.5 < z < 1 

1< z < 2 

2 < z < 5 



212 
43 
70 
68 
15 



74 
26 
24 
14 
3 











Aird et al. (2010 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



223 
44 
97 
68 
4 



99 
32 
45 
16 




16 
3 
8 
5 




16 

8 
7 
2 












DB10 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



230 
74 
79 
67 
10 



91 
46 
30 
14 
1 



19 
3 
6 
9 
1 



14 

6 

6 

3 













Treister et al. (2009 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



191 
51 
65 
65 
11 



77 
35 
26 
14 
1 



4 
2 
1 
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Model 



Table 2 — Continued 



Redshift Range 



All AGNs 



Compton Thick AGNs 



6-10 keV 10-30 keV 30-60 keV 6-10 keV 10-30 keV 30-60 keV 



AEGIS (0.5 deg 2 ) 



Uedaet al. (2003 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



171 
34 
57 
61 
11 



62 
22 
22 
13 
2 











La Franca et al. (2005 



< z < 5 
< z < 0.5 
0.5 < z < 1 

1< z < 2 

2 < z < 5 



156 
26 
51 
58 
12 



59 
16 
20 
15 
3 











Aird et al. (2010 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



144 
24 
58 
52 
5 



71 
18 
32 
17 




17 

3 
7 
7 




15 
5 
7 
3 












DB10 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



184 
43 
60 
68 
13 



71 
28 
25 
16 
2 



19 
2 
4 
10 
2 



11 

3 
4 
4 












Treister et al. (2009 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



151 
29 
49 
60 
13 



59 
21 
22 
15 
1 
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Model 



Table 2 — Continued 



Redshift Range 



All AGNs 



Compton Thick AGNs 



6-10 keV 10-30 keV 30-60 keV 6-10 keV 10-30 keV 30-60 keV 



ECDF-S (0.25 deg 2 ) 



Uedaet al. (2003 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



145 
26 
48 
54 
12 



54 
17 
20 
13 
2 











La Franca et al. (2005 



< z < 5 
< z < 0.5 
0.5 < z < 1 

1< z < 2 

2 < z < 5 



129 
20 

42 
49 
11 



51 
13 
18 
15 
2 











Aird et al. (2010 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



111 

17 

43 
43 
5 



58 
14 
26 
16 
1 



16 
2 
6 
7 
1 



13 

4 
6 
3 












DB10 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



156 

32 

50 

61 

14 



62 
21 
22 
17 
2 



17 
1 

4 
10 
3 



10 
3 
3 
4 












Treister et al. (2009 



< z < 5 
< z < 0.5 
0.5 < z < 1 
1< z < 2 
2 < z < 5 



125 
21 
40 
51 
13 



52 
16 
19 
15 

2 











Table 2 — Continued 



Model Redshift Range All AGNs Compton Thick AGNs 

6-10 keV 10-30 keV 30-60 keV 6-10 keV 10-30 keV 30-60 keV 



GOODS (0.089 deg 2 ) 



Ueda et al. C2003) 


< z < 5 


106 


45 


1 


6 


4 







< z < 0.5 


17 


11 


1 


1 


2 







0.5 < z < 1 


34 


16 





2 


2 







1< z < 2 


41 


14 





3 


1 







2 < z < 5 


11 


2 


Q 


1 








La Franca et al. (20051 


< z < 5 


94 


41 


1 


5 


4 







< z < 0.5 


13 


9 


1 


1 


1 







0.5 < z < 1 


31 


14 





1 


1 







1< z < 2 


36 


13 





3 


1 







2 < z < 5 


9 


2 





1 








Aird et al. (2010) 


< z < 5 


75 


41 


1 


13 


10 







< z < 0.5 


10 


9 


1 


2 


3 







0.5 < z < 1 


27 


17 





5 


5 







1< z < 2 


30 


13 





6 


3 







2 < z < 5 


5 


1 





1 








DB10 


< z < 5 


115 


51 


1 


15 


8 







< z < 0.5 


21 


14 


1 


1 


2 







0.5 < z < 1 


35 


17 





2 


2 







1< z < 2 


46 


17 





8 


4 







2 < z < 5 


13 


2 





3 








Treister et al. (2009) 


< z < 5 


87 


42 


1 


2 


2 







< z < 0.5 


13 


10 


1 





1 







0.5 < z < 1 


28 


15 








1 







1< z < 2 


36 


14 





1 


1 







2 < z < 5 


10 


2 















Note. — These are 4cr detections for the half shift surveys. Numbers for other tilings/surveys are available from the 
authors. The models are described in Sect. 12.41 The three fixed Compton thick fraction models utilize the HXLFs of 
Ueda et al. (2003; / CT = 0.3), La Franca et al. (2005; / CT = 0.4) and Aird et al. (2010; / CT = 0.5). The DB10 model 
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is the 



Draper fc Ballantvnd (|2010T ) model of Compton thick evolution and makes use of the lUeda et al. ( 2003 ) HXLF . 



Finally, the the most conservative v iew of the Compton thick population is based on the model of Treister et al. ( 2009|) 
and also uses the Ueda et al. I d2003h HXLF. Numbers have been rounded, so the sum of the individual redshift ranges 
may not always equal the value in the < z < 5 row. 



